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CONTRIBUTION TO THE THEORY OF STATIONARY SEPARATION AREAS 


G. I. Taganov 

Experience shown that the pattern of steady flow around /3 t 
a bluff body with a closed separation area is disrupted at 
subsonic current velocity at R numbers from 10 1 to 10 2 when the 
flow is laminar in nature. In addition, experience demonstrates 
that, at supersonic current speed, there is a stable, steady 
flow pattern if there are laminar stagnation zones adjacent to 
the body (a stagnation zone behind a reverse ledge on the body's 
surface, a stagnation zone in front of the smooth ledge on the 
body’s surface, a forward separation area formed by the tip of 
the centerbody, a stagnation zone formed when there is a 
descending junp on the body’s surface) at high R numbers on the 
order of 10 4 -10 6 . 

Thus, experience indicates that, for certain ranges of change 
in M and R numbers under specific boundary conditions, steady 
solutions to equations for the movement of n viscous liquid exist 
and are reliable. Outside these ranges and under different 
boundary conditions, flow around a body with a closed separation 
area is more (Earcon’s vortex street at M<<1) or less (pulsing in 
close wake behind a body at M>1) clearly definable as steady in 
nature. There is still no theoretical proof of the existence of 
stable steady flows with separation areas in these ranges. 

Here we attempt to find the region of existence of possible 
steady flows with a closed separation area in a range of R 
numbers such that flow in a viscous nixing area can he described 
by Prnndtl’s equations. Limit conditions for flow within the 
separation area are set so that the flow pattern within the area 


^Numbers in the margin indicate pagination in the foreign text. 
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is considerably simplified and it becomes possible to use 
hydrodynamic analytical methods. 

The first part of this article (points 1-4) studies the field 
of possible steady flows for an incompressible fluid. It is 
shown that only under special boundary conditions within the 
separation area (ideal dissipator) does flow near a smooth slab 
approach Kirchhoff’s flow with quiescent fluid within the zone if 
R approaches infinity. The drag coefficient for the "slab + 
ideal dissipator" system is c (n + 4 ) , i.e. nearly half the 

value obtained by Kirchhoff for an ideal fluid. 

Qualitative study of the field of possible steady flows in 
plane c^R has made it possible to detect the existence of a 
certain region which has degenerated into a line with an upper 
limit on the R number of about 10 1 2 . In this region steady flows 
possess a particular flow configuration if vortical attachment is 
not viscous. 

The existence oi a connection between flow configuration in 
the region of non-viscous vortical attachment and the reliability 
of the steady solution is traced in part two (points 6-7), both 
for individual solutions obtained by methods of the linear theory 
of hydrodynamic stability and for well-known experimental 
material obtained over a wide range of R numbers at sub- and 
supersonic current velocities. This review makes it possible to 
formulate a rule for selecting stable steady flows with 
separation areas and to use it to analyze comparable flows — 
both laminar and,, in some cases, turbulent. 

1. Degenerate flow in a separation zone. We will consider 

the pattern of a steady flow with closed separation urea behind a 
body past which a current of incompressible fluid flows at R 
numbers such os R=u — d/v (u^ is the velocity of the unperturbed 
current, d - the transverse dimension of the body), at which the 




thickness of the boundary layer is less than the size of the body 
(cf. figure 1). 

We can identify the following typical lines and areas in the 
upper half-plane of the flow* 

The dividing streamline of the current with fluid segments AB 
and BC and a segment of the body boundary AC; regions: 

1 — outer potential current; 2 — outer viscous boundary layer 
and wake; 3 — inner viscous mixing layer; 4 — viscous flow with 
nearly constant vorticity; 5 -- viscous boundary layer near the 
rigid boundary of the body v/ith a vorticity opposite in sign; 

6 — attachment* At points on segment AB of the dividing 
streamline the fluid has a velocity of u^, and friction 3 tress 
is applied clockwise to the fluid moving within the zone. In 
segment CB, fluid velocity is u 0 ; friction stress equals zero. 

In segment CA, velocity equals zero, and friction stress applied 
to the fluid acts counterclockwise. 


During 
segment AB 
the closed 


s t eady 
should 
area . 


flow in a separation area, friction forces in 
equal the viscous dissipation of energy within 


' A» 




( 1 . 1 ) 


Here dl is an arc element, dj is an element 
section of the area, F is the area of the cross 
by line ABCA. 


of the crocs 
section limited 



Analysis of values introduced into the integrated expression 
on the right side of equation (1-1) shows that, if R numbers 
increase, when the thickness of regions 3 and 5 decreases, 
velocity gradients in then noticeably exceed the velocity 
gradient in the core of the flow. Most of the dissipation occurs 
in these regions, and most energy dissipates in that part of 
region 3 which is adjacent to section AB of the streanline. 


He will try to find conditions under which the pattern of 
flow within a separation area may be fundamentally simplified. 
To obtain such a degenerate flow within the separation area, we 
need to: 


1) prevent energy dissipation in region 4; 

2) minimize dissipation in region 5; 

(3 ) 

3) precisely determine energy dissipation (E that 

part of region 3 which is adjacent to the st>eamline AB , 

4) concentrate dissipation 


• *L 




AB 


■ Fleur* 2 


inadequate to provide a balance in segment CB of region 3, by 
causing dissipation in an additional element — an energy 
dissipator. It is assumed that this additional element, firmly 
connected to the body, creates a region with high velocity 
gradients within the separation area, in which dissipation 
inadequate for balance (4) occurs and that it accepts and 
transmits to the body the momentum of the backward jet moving 
toward the body from the attachment region. The shape and 
dimensions of the dissipator are immaterial to further 
discussion — it is considered to be a black box with the ability 
to reduce the value of Bernoulli constants cor responding to 
different backward jet streamlines to assigned values. 
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Figure 2 shows a flow arising in the presence of an ideal 
dissipator, which dissipates almost all the kinetic energy of the 
backward jot and assumes the entire momentum of the backward jet 
gained in segment AB of the dividing streamline. 

It is easy to see that, in this case, degenerate flow ZJL 

(without a circulation core) occurs at constant pressure over the 
entire region limited by the dividing streamline, except for the 
attachment region. 

A viscous mixing layer develops here at constant pressure and 
constant outer flow velocity. The tangential component of inner 
flow velocity on the border of the viscous nixing layer equals 
zero. Note that there are precise self-similar solutions under 
conditions close to those considered [1, 2]. We will calculate 
the residual kinetic energy contained in the backward jet, which 
must be dissipated in the dissipator, as well as the backward jet 
momentum which the dissipator must assume. 

Let us calculate parameters of the inner part of the viscous 
nixing layer in cross section S, located at the end of the 
separation area, but in front of the attachment region in wnich 
the flow will be regarded as nonvisr.ous and eddying. 1 

The force of the outer current at the dividing streamline 

(1.1) is partially dissipated when the quiescent fluid 

accelerates and is partially used up in creating motion whose 

kinetic energy in cross section S equals c 

VsP ] tt l {n)dn 


i It is possible to regard flow in the stagnation point region 
as nonviscous and vortical because of analysis of Navier— 
Stokes equations in this area. (This has been shown, for 
example, in a paper by V. V. Sychev at the 8th Symposium on 
Modern Problems of Fluid and Ga3 Mechanics. Poland, Tarda, 
18-23 September, 1967.) 
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where n is the normal to the dividing streamline. This residual 
kinetic energy should dissipate in the dissipater with an error 
on the order of (v*/u p ) 2 , where v* is the velocity of the fluid 
drawn into the viscous nixing layer from the quiescent fluid in 
the separation area and the level of dissipation in the section 
between cross section S and the dissipator. 

Let the non-self-s imi lar viscous nixing layer in cross- 
section 1 have a certain momentum. It will remain the same in 
cross-sect i on 2 dowstrean as well, but will increase in the inner 
part of the viscous mixing layer in cross-section 2 and decrease 
by the same amount in the outer part, since the momentum gained 
by the inner part of the mixing layer at the area of dividing 
streamline ^-l^ be 

P (l)dl 

i, 

and the momentum lost by the outer current will be AV. 


If the thickness of the mome 
separation area in the outer cur 
cross section S, the momentum ga 
current will equal 


n turn 

1 O £ 

; t at 

the beginning of the 

rent 

i s 

othe 

r than zero, then in 

ined 

by 

the 

inner part of the 


0 CD «J 

[p 5 uyn ]_=[p$ »(«k — »)**],— [pj »(«** — B) At] o 

—eo 0 0 


( 1 . 2 ) 


where is the velocity at the outer boundary of the nixing 
layer. If we use the concept of momentum loss thickness and 

introduce the concept of the momentum gain thickness of the inner 
part of the viscous mixing layer + determined by the /6 

relationship 

c 

r= p J urdn (1.3) 
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then from (1.2) we have 


<6 + °*)s = (6 s *)s- (6**)o 


(1.4) 


The combination of the constant pressure region behind the 
body and the inner potential flow region was studied using the 
theory of cavitation flows. At present, there is a mathemati- 
cally noncontradictory pattern for flow with a backward jet 
(Efros’ flow [3]), which ensures combination of the constant 
pressure zone 



Figure 3 


1 


behind a body and inner potential flow over a wide range of 
static pressures within the zone (cf. figure 3). However, a 
hydrodynanically noncontradictory outer potential flow around the 
body with constant pressure behind it must be plotted, in which 
viscous effects within the closed separation area and viscous 
effects in the outer part of the nixing layer (displacement 
effect) would be accounted for. 


In approximating momentum equations of the boundary layer 
theory, flow in the inner part of the viscous mixing layer with 
degenerate flow within the area is equivalent to flow in a fluid 
layer with Bernoulli constant for outer potential flow and with a 
thickness equal to nomenturj gain thickness d + *«. 

Consequently, the thickness of a backward jet with a velocity 
constant in terms of cross section and corresponding to the 
pressure in the constant pressure zone in ideal fluid model <5° 
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must equal the momentum assumption thickness in cross-section S 
of the viscous flow: 


6 °=( S + ")6 


( 1 . 5 ) 


Tractive force T, acting on on ideal dissipator in a viscous 
degenerate flow 2 , must equal tractive force T° , which acts on 
flow within the constant pressure zone in on ideal fluid model of 
flow with a backward jet. 

In viscous degenerate flow, the ideal dissipator is a flow of 
impulse, but not a flow of mass. In the ideal fluid potential 
flow which models the subject viscous flow, sink flow of impulse 
which equals the power of the sink flow of impulse in an ideal 
dissipator, is, of course, sink flow of mass. 

Presure drag X t acting on a body in viscous degenerate flow 
must equal pressure drag acting on a body in on ideal fluid model 
X 0 

A i • 

lhe displacing effect a viscous boundary 1 ayer developed at a 
rigid profile has on potential flow is, as we know, determined by 
the arrangement on the surface of the profile of hypothetical 
sources, which ensure that potential flow deviates by 
displacement thickness 5V. Hypothetical sources in an actual 
viscous flow do not have a current of impulse. In on ideal fluid 
model, the current of mass from these sources also bears the 
current of impulse. 


In the case considered — flow with a free boundary — /I 

current from sources arranged close to the dividing streamline 
deforms not only exterior potential flow, but also the contour of 
the constant pressure zone. To obtain from Efros’ flow an ideal 


2 Here and henceforth, T, Xj , and X represent total forces 

valid for two flow regions divided by the plane of symmetry. 
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fluid model in which viscous effects in the internal and external 

part of the viscous layer are accounted for in case of degenerate 

flow in the separation area, sources must be placed close to the 

free boundary between cross sections 0 and S. These sources must 

ensure that a random cross section 1 has a mass flow rate of 

P u, (d ** + 6 *), where 6 and C # are, respectively, momentum 
k » e t e 

gain thickness and displacement thickness of the outer part of 
the viscous mixing layer in this cross-section (cf. figure 4*) 

Figure 4 

C>nr. 4 

Indeed, let us extend line L from point N, the dividing 
streamline of Efros* current (cf. figure 3) in cross section S, 
to point 0, so that the current of impulse in each section of 
each layer limited by line L and the free boundary will equal 
PV (* + **) 1 an ^’ correspondingly, the current of moss will equal 
pu k (5 + ^) 1 . 

Then for the external flow, line L will become the line of 
sink flows through which nassflows equal to {-pu ^ can 
successfully pass. To account for the displacing effect of the 
external part of the viscous layer, we must arrange along this 
line both sources which compensate for sinkflow on the line with 
a flow rote of toward cross section 1, and sources 

which account for the displacing effect of the external part of 
the viscous mixing layer with a masoflow of toward 

cross section 1, i.e. with a total massflow toward section 1 
equal to pu^((J + ^*^ + <5^)^. 

If in the first approximation we disregard the effect of 
viscous layer displacement on the configuration of the constant 
pressure zone, then error in calculating transverse dimension b 
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?'A 



of the constant pressure zone will be on the order of 

6 */b R. ~ 1/Z . 
e b 

In direct use of Efros’ flow in a first approximation of 
external potential flow near ■ body with degenerate flow in the 
separation area, as the preceding indicates, the outer part of 
the viscous nixing layer is characterized by negative 
displacement thickness, equal to the momentum gain thickness of 
the inner flow, In other words, an additional error in 

determining the transverse dimension of the constant pressure 
zone, on the order of about R b _1/2 , is permitted. 

We will now determine in the first approximation (without 
regard for actual displacement) c^ for a "flat slab-ideal 
dissipator system as a function of H number, assuming that the 
initial momentum loss thickness (d**) 0 =0 and flow in the viscous 
nixing layer is self-similar. 

In this case, with regard for (1.4), tractive force acting on 
the ideal dissipator is by definition equal to the total backward 
jet momentum (since the axial component of the velocity of the 
fluid from the dissipator equals zero) 


%T = p^(6* K )6 


( 1 . 6 ) 


Pressure drag acting on a bluff body, is about equal to 
X 1£ . Pressure drag acting on a body in Efros’ flow is related to 
backward jet thickness by a ratio obtained from the equation 
for momentum (cf. [4]) 



(1.7) 


With regard for tractive force acting on the flow of fluid /J3 
in a backward jet l/2T E =pu^d E , the coefficient for total force of 
drag acting on the "body + backward jet" system will be 
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(1-8) 


rarMst-s-xsc 


2X 

pu«,y 


-«?(£> 


and the coefficient of pressure drag acting on the body is, from 
(1.7) 

*>Fj 4 , «M tt * 
u« 


, O’E/ A \ «i 

Cxi = 4 — i 1 H ) — 

a \ u« / u a 

where d is the size of the body. 


(l.P) 


The dimensionless velocity at the current dividing line in a 
self-similar solution [1, 2] will be 

U = — = 0.5L7 


and friction stress on the current dividing line (from [2]) will 
be 


ip = 0.2y’{iP“* 3 / l 


( 1 . 10 ) 


where 1 is the length of the arc along the dividing s t reanl ine 
fron section 0 to the subject section. 


The configuration of the flow, and then also for a given R 
number, can be determined from (1.5), since the length of cavity 
lj t as a func* ’on - of backward jet thickness is known from the 
solution of the problem of streamlining a body with backward jet 
fj ( 1 ) • From the self— similar solution for a viscous mixing 
layer, we can obtain the function (<S + **)„=(<J**) s = ( 1 Ic ) For a 

given R nunber. 


Integrating (1.10) along the dividing streamline to cross 
section S, located distance 1^ from the beginning of the nixing 
layer, we obtain an equation for the momentum loss thickness of 
the external part of the viscous mixing layer 


pB k *(6**) B = 0/iViipur^ 


( 1 . 11 ) 



Substituting the equation for from (1.11) into (1.8), 
given (1.5) and (1.6), we find 

c* « 1.6 h' ; 'U k '''R-% Ik = Ik/ d, u k = u k /u x (1.12) 

1^ is the dimensionless length of the cavity in an Efros 
is the dimensionless velocity nt the free boundary. 

From (1.8) and (1.9) we have c ~c ./(1+u,). As we know [4], 

.. X X X K 

when lj c approaches infinity, then in an Efros flow near a flat 

plate 2ffg/d^TT/s (tt+ 4) , and c . approaches 2 tt/(tt+ 4), about 0.88, 

i.e. it tends toward a Kirchoff value for drag coefficient. It 

is easy to see that (1.5) is satisfied if R approaches infinity. 

Consequently, flow in a separation area with an ideal dissipator 

approaches a Kirchhoff flow when R approaches infinity, and, 

according to (1.8), when u. approaches u , the drag coefficient 

K ** 

for the "plote + ideal dissipator” system, approaches c^ = , n/(TT + 4) 
or about 0.44, i.e. a value half that obtained by Kirchoff for a 
separation area with a quiescent fluid at w=0. This is an 
accurate result, since error from disregarding the effect of 
displacement approaches zero if the increase in the size of the 
constant pressure zone is unlimited. 


Here 

flow, 


The limitation on 
flows in terms of R n 
PQp is total pressure 
vortical flow in the 
used in the region of 
streamline and to the 
follows that 


the region of existence 
umber imposes the condit 
at the dividing streaml 
attachment region. From 
nonviscous attachment t 
streamline for the inne 
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(1 

. 13) 


When R approaches infinity and U=0.587, Pq^=0.345. As R 
decreases, rises, and for a certain limit R, value R^^, 
becomes zero. The value u^» corresponding to this limit R 
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number is found from (1.13) when Pq^ = 0. For U = 0.587, 1*23 

and the corresponding value for pressure resistance in the 


separation area is 




= 1 



Let us calculate R, . and c corresponding to the limit 

1 IQ X 

steady flow around a slab with degenerate flow in a separation 
area. Using results calculated to solve the problem of 
streamlining a slab according to Efros* system [4], when p^= 
-0.52, we have 1^ equal to about 15, equal to about 1.23, c^ 
equal to about 0.03+0.52 and, from (1.8) and (1.9), equal to 
about 0.62. 


Then, substituting these numerical values into (1.12), we 

obtain an R. . of about 120. Thus, in plane c ,R (cf. figure 5) 
1 in * 

all steady flows around a slab with degenerate flow in the 
separation area lie on curve DC, extending from un about 

10 2 to R-+~. Along this curve, c v for the "slab + ideal 
diss'pator" system varies from about 0.62 to about 0.44. 


2. Nondegeneruto circulation flow in a separation area. Let 
us now consider q flow in a separation zone with an imperfect 
dissipator which permits the existence of circulation motion in a 
separation area. Batchelor [5] showed that a flat flow of u 
nonviscous fluid with constant vortex value over the entire 
circulating flow region, except for a thin layer close to the 
boundary of the region whose thickness approaches zero when 
R+~ . 
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Let 

constant 

velocity 

external 

develops 

dividing 

boundary 


us consider a boundary separating circulation flow with 
vortex from external flow with a zero vortex value. If 
on the boundary line is less than velocity of the 
potential flow at this line, a viscous nixing layer 
if p does not equal zero. Friction stress on the 
streamline of the viscous nixing layer causes the 
layer of the circulation flow to accelerate. 


Let us place a dissipator across the streamlines of this 
circulation flow boundary layer which is being accelerated and, 
consequently, is experiencing a positive increase in Bernoulli 
constant. Let the dissipator take not all, but only the 
additional momentum gained by this boundary layer when exposecj to 
friction stresses on the dividing streamline and let it release 
fluid with the same Bernoulli constant distribution across the 
streamlines os was in this layer before it accelerate near the 
dividing streamline. The presence of a dissipator with these 
properties in the separation area with circulation notion /IQ 

makes it possible to satisfy equations for energy and momentum. 

The velocity of the external potential flow u. and the 

k 

velocity at the boundary layer of the circulation flow with 
constant vortex u^-u^/u^ along the common boundary are related by 
the equation 


A = a** — u 2 = const 


( 2 . 1 ) 


resulting from the equality of 
and internal flow at each point 
constant introduced into (2.1) 


s ta 
on 


tic pressure in the external 
the boundary [6]. The 
be used as an independent 


can 



it 
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;j 

parameter defining flow and, accordingly, dissipator properties, 3 ! 



An ideal fluid model describing the external potential flow 
near a body with a separation area with circulation flow in it 
and accounting for all viscous effects is built similar to that 
in point 1. The dissipator in a viscous flow will be a sink flow 
of mass; in the ideal fluid model, a sink flow of impulse with 
something like the power of a sink flow of mass will correspond 
to it. Strictly speaking, the location of this sink flow should 
correspond to the location of the dissipator in the separation 
zone (in this it differs from degenerate flow). 

Figure 6 


The currents of mass in any radial crocs section of 
circulation flow in an ideal fluid model are identical (cf. 
figure G) . Calculation of parameters for flow in a viscous 
mixing layer (defining U(l), tffc*(l), <5 e ^ ( 1 ) ) is complicated 
because the difference in velocities on the boundaries of the 
layer, as well as static pressure, varies along the viscous 
layer. 

There are also calculating problems in finding flow in an 
ideal fluid model, resulting because one region of flow satisfies 
Laplace’s equation, while the other -- the inner — satisfies 

3 If a portion of the friction surface of the body acts as a 
dissipator, then parameter A is defined by the condition 

valid for a boundary layer periodic in terns of 1 at high R 
numbers. 




17 



4, + f , U *%* j,T<« iW W» 




« *r^+r: ^ r— : * »vs — *■ ar ■ 


Poisson's equation. We con approximately replace flow with 
constant vortex equivalent to it in terns of current of impulse 
by M. A. Luvrent ’ yev ’ s potential ring [7] and seek a solution to 
the problem in the class of harmonic functions. Numerical 
calculations are carried out in this approxinat ion. 


r-r> ■*>-*-*■* £/ J-. J 

0 


a ‘ 

ll 

l 

:1 


3. The Region of stationary flow in plane c^R. Without data 
fron numerical calculations on flow with circulation motion 
within a separation area, we will try to qualitatively evaluate 
the location in plane c^R of stationary flows which satisfy 
various values for parameter A (cf. figure 5). Steady flows with 
an ideal dissipator in the separation area for which u^eO 
correspond to points on curve BC. According to calculation, 
parameter A varies monotonical ly along the curve, according to 
calculation, from A 0 of about 1.52 to A c equal to 1. The lower 
boundary of the region of steady flow with circulation motion in 
the separation zone in plane c^.R must be curve AD, which 
corresponds to a maximally weak dissipator (friction on the back 
side of the slab). It is easy to show that the drag of a slab 
streamlined in the separation area corresponding to points on 
curve AD is on the order of the friction drag in the seme slab, 
post which a current flows at the same R number, but below zero 
incidence, i.e. curve AD is close to curve L, described by the 
equation c =2( 1 . 3'28R_ /2 ) . 0 

fp l u(o — Erjrfr]* 0 

— <x 

Parameter A on curve AD cannot become zero when R number is 
finite, since the finite friction impulse on the back side / j x 

of the slab which decelerates the circulation flow within the 
area and the absence in this case (u^-u^O) of cn accelerating 
impulse from the internal flow at the dividing streamline would 
violate the equality •Uj { 2 -u b 2 =0 and create o new circulation flow 
with a A> 0 . Parameter A can become zero only if the 
decelerating friction momentum on the back side of the slab 
becomes zero, i.e. if R+- (point D). If R+-, A, „hich describes 
dissipator efficiency, passes through all values between zero at 
point D(c x +0) and unity at point C ( c *u/u+4 ) . Keeping this in 
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nind and relying on the distribution of A calculated along curve 

DC, we can expect that the behavior of lines A = const in plane 

c R io close to that presented in figure 5. 
x 

We will now determine how the shape of the separation area 
close to the slab and the flow close to it vary when R-~ as 
dissipator effiency decreases, i.e. during movement from point C 
toward point D. Point C represents an infinitely long separation 
area, as indicated above. The unlimited size of the separation 
area for points located between C and D when R-*~ results from the 
finite nature of at these points. 

Even an approach to point D along curve AD does not 
contradict the assertion that a flow corresponding to point D 
also belongs to a class of flows with unlimited separation area 
dimensions. The length of the separation area must increase in 
flow corresponding to points on curve AD if R+~; otherwise the 
accelerating momentum, applied to the circulation flow_at the 
dividing streamline approaching zero as R t n j t a "' n b 2)1 lt* would 

decrease more quickly than the impulse of ^friction on the back 
side of the slob approaching zero when R 

Flows with. a vorticity 0 of about u w /b, where b is the 
transverse dimension of the area, correspond to all points of 
segment CD, except degenerate flow in the area corresponding to 
point C. Since the size of the area increases without limit when 
r->-<*> , flow within the separation zone at a distance on the order 
of d from the slab approaches local potential. 


Consequently, at the tangential discontinuity separating the 
separation area, it is not vorticity, but the Bernoulli constant 
that is disrupted if A>0. 


The shape of the separation area in plane xy when the fluid 
in the separation zone is quiescent, taken from Sirchoff*3 well 
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known solution, is given in figure 7. 

I' 

/ 

Figure 7 



Oar. 7 




If we assume the existence in this contour of a potential 
flow with A>0, then pressure at the stagnation point on the back 
side of the slab will become greater than p^. It will also 
increase m the angular region of the subject contour adjacent to 
the slob. To fulfill the condition of equality of static 
pressures for both sides of the contour, it must deform — the 
curvature of the contour in the angular region must decrease as 
Compared with initial curvature. As parameter A decreases, the 
curvature of the contour should diminish (cf. figure 7). jfe can 
find the shape of the separation area contour from the solution 
of the problem of streamlining a slab with two oncoming clean 
currents for which the difference in Bernoulli constants is set 
when static pressures are equal and the fluid contour is 
impenetrable. 

» 

Streamlining a slab with two oncoming clean currents with 
identical Bernoulli constants (A=0), which provides a local flow 
pattern near the slab, a pattern corresponding to point D on 
plone c^.R, is worth discussing. 

It is easy to be sure that the conditions of the problem are 
satisfied by a potential flow whose flow function is 

—i! (3-D 

' °V dz ) 

Coefficient k from the condition that velocity is limited at 
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an infinite (to the scale of dimension d of the slab) distance 
from the slab should approach zero, and flow near the slab should 
become quiescent (of. figure 8, a). 

Thus the limit condition for the flow of a viscous fluid /12 
near a slab with a separation zone when R-~ is in fact, as 
Batchelor first stated [6], by c^O. However, the flow pattern 
turned out to be quite different from that assumed. Nonetheless, 
the resulting flow near the slab satisfies hydrodynamic equations 
and provides another description of Dalamber's paradox for flow 
around a slab: it accounts for the need to form a separation 

area because of the effect of viscosity, which is usually 
disregarded (cf. figure 8, b). 

Let us turn now to the other end of region ABCD — to small R 
numbers . 


Curve AB is the boundary of the region of possible steady 
flews. Its points represent limit flows for which the equation 

p 0p“ (p 0p /l/2pu^ 2 in the region of nonviscous vortex 
attachment becomes nero, 


Using calculations of the change in p Qp along curve BC, which 
varies monotonical ly in R number from (F ) B =0 to (p )c=0.345, 
and taking into account that p Qp approaches 1 when parameter 
epproaches zero, we can construct the preferred pattern for 
locating lines p Qp =' const in the region ABCD (cf. figure 5). 

Of course, when analyzing flows in the region of low R 
numbers, we must keep in mind that the original assumption that 
flow in a viscous nixing layer con be described by Prandtl’s 
equation and that circulation flow within an area can be 
described by Eiler’s equations ore invalidated for a certain 

value R nin , and if R<R rain , flow is described by Novier-S tokes 
equations. 
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4. Unique features of flow in the region of nonviccouc 

vortical attachment when p Qp =0. Figure 9 presents possible 

alternatives for pressure distribution 

5 = £ P°° 

l h pu« 2 

along the plane of f lov; symmetry for nonviscous attachment of 
vortical currents in a separation area behind a body, which 
correspond to various points in region A BCD and its boundaries. 
At the maxima for these curves, we have P~Pgp> and the position 
of the maximum along the x-axis determines the position of the 
stagnation point behind the separation area along the y-axis. 



JK. 


Figure 8 


If, for flows 1 and 2, for which p 0o >0, the stagnation point 
lies farthest from the separation zone, then for flow 3, for 
which PQp“°> three types of pressure distribution curves along 
the plane of symmetry are conceivable: 3a — pressure maximum 

lies within the flow, but there must be at least one pressure 
minimum downstream of it; 3b — pressure maximum is achieved 
within the flow, but downstream of it, pressure is constant and 
equals pressure at an infinite distance; 3c — the stagnation 
point lies an infinite distance from the stagnation zone. Using 
A. A. Nikol’ skiy ' s theorem [8] on the nonotonic nature of the 
change in the slope of the velocity vector during movement of a 
gas along the equal pressure line in a flat subsonic vortical 
flow, we can show that flow conditions described by curves 3a and 
3b are impossible in the cose we are discussing. 
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If we extend the constant pressure line near the assumed j 

minimum pressure lying downstream of the pressure maximum j 

corresponding to condition 3a, then the ends of this constant j 

pressure curve lying on the plane of symmetry will be / 1 3 

represented by identical velocity vector slope values, and this 
will contradict Nikol’skiy’s theorem, which requires that 
velocity vector slope increase if the length of tne equal 
pressure line is finite. 


If we extend the constant pressure line from the point 
corresponding to the assumed pressure maximum in 3b, then it goes 
into infinity, where the velocity vector slope is also equal to 
the velocity vector slope at the plane of symmetry, i.e. we will 
arrive at the contradiction already discussed. 

Consequently, when ^^=0, the nonviscous vortical flow in the 
attachment region has a configuration (cf. figure 10) such that 
the stagnation point does not lie in the finite region of the 
flow. In this flow configuration, the dividing streamline 
(denoted in figure 10 by the thick line) limits the thin tongue 
of the nonviscous backflow extending downstream to infinity. 



Thus, points on boundary AB of region ABCD (cf. figure 5) 
represent steady stalled flows whose stagnation point does not 
lie in the finite region of the flow when attachment is 



/ 
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nonviscous and vortical. In other words, these flows have a 
configuration in attachment regions significantly different from 
that of a flow with a stagnation point in the finite attachment 
region, which is typical of flows filling the entire region AGCD. 

The existence of a connection between the configuration of a 
flow in the nonviscous attachment region and the stability of a 
steady flow with a separation area is traced in points 6 and 7. 

5. Steady flows with separations areas when H>1. When M>1, 
only part of the total loss in momentum caused by base drag lies 
in the vortex wake. Most of the momentum los 3 is caused by the 
wove drag of the fluid contour of separation area X 2 and lies in 
waves moving downstream from the separation area boundary. 

Base drag of a flat body of finite length will be 

'/»*■= + 

The coefficient of bottom drag related to the velocity head 
of an unperturbed flow and to the finite dimension d of the body 
in the base shear region will then be 

«=24.**/d4-c„ 

If ( d ££ ) q = 0 and flow in woke behind 
considered nonviscous and vortical, then 
contour ofthe separation zone can be wri 

C„ = e,„ (1— 25„*/<0 

Here c j{2f is the coefficient for the wave drag of the fluid 
contour of the separation area related to the height of the fluid 
ledge d-2ff^*. 

clearly, the wave portion of base drag depends on the shape 
of the separation zone’s fluid contour and wove displacement 
thickness. The development of wave drag, which is a function of 
the separation zone’s fluid contour, and separation zone contour 


the shear region is 
wave drag in the fluid 
tten as 
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as a function of circulation flow intensity within the area make 
it much more complicated to find possible steady flows in plane 
c^R over a wide range of changes in dissipator properties* 
However, for degenerate flow within the area, i.e. if an ideal 
dissipator is present, the problem is considerably simpler and 
possible constant-property flows can be found. We will limit 
ourselves here to discussing the case of a degenerate flow within 
the area when and 0. The self-similar precise solution 

for a laminar layer on the dividng boundary nixing a jet of 
compressed gas with a quiescent gas when yp = const from [2, 14] 

gives us 



Uk/C * 



(5.1) 


Total pressure on the dividing strearaling in the non- / 14 

viscous attachment region is calculated from 

■ — = — . *>(*)«[ i— H-ivT’ Y = — 

Pm p(Mp(M L Y+i J (5.2) 


Static pressure within a separation area which has a straight 
boundary is constant over the entire area; d+0 and d^fc/d^O 

when R-h*, while c x Hc x 2f* Static pressure in the area, which 
determines base drag is in this case obtained from the ratio 


<) 


•ltd 
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Figure 11 


(5.3) 


From (5.1)-(5.3) for a given we can calculate p^/p^ which 

corresponds to possible steady flows when R-*~ in plane c x R (cf. 
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figure 11). Point D represents p Q /p..>l, point C corresponds to j 

the limit flow with a separation area P 0 p /p„=l and limit values j 

(p, /p ) , . with c , . . i 

i 

* 

We can show that flow configuration in a nonviscous vortical j 

attachment region with infinitely long return flow tongue \ 

corresponds only to limit flow PQp/p )W = l» just as to limit flow 
p Q p=0 when M+0. All remaining steady flows with POp ^ 1 contain ! 

a stagnation point in the finite region of the nonviscous 
vortical attachment region. 

Using the theorem of monotonic change in velocity vector 
slope along a constant pressure line located in the subsonic 
region of a flow adjacent to the plane of symmetry, i.e. assuming 
that we can always select a constant pressure line close enough 
to the proposed minimum pressure that this line is within the 
subsonic region of the flow, we will establish that such a 
pressure minimum is impossible, i.e. that condition 3a (cf. 
figure 9) is impossible. 

The hypothesis that condition 3b can exist implies that 
pressure must be constant over the entire subsonic region of a 
flow located downstream of a constant pressure line extended from 
the point of the assumed pressure maximum. But this in turn 
requires that both pressure in the supersonic region and velocity 
vector slope be constant. During flow around a flat 

M~f 

Figure 12 


body of infinite length with curved shock waves upstream, this 
requirement cannot generally be fulfilled. The configuration of 
flow in a nonviscous vortical attachment region, corresponding to 
limit flow with P 0 p /P«, = 1 > is shown in figure 12 . 
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6. On the existence of a relationship between configuration j 

of flow in a nonviscous vortical attachment region and flow j 

stability. As we know, study of the stability of a given steady j 

flow in terns of low velocity and pressure perturbations } 

depending on tine and coordinates presents a difficult j 

hydrodynamic problem which, within the framework of the linear \ 

theory, can be solved only for individual simple configurations \ 

of viscous steady flows. Study of the rtability of viscous j 

steady flow with a separation area of complex configuration also f 

presents problems which are yet to be overcome. However, study j 

of the hydrodynamic stability of nonviscous steady fluid flows 
involves fewer mathematical difficulties. Therefore, we do have 
some results, obtained with the theory of hydrodynamic stability 
for several nonviscous flows of incompressible fluids of complex 
configuration with partially free and partially rigid boundaries 
which contain regions of jet attachment or connection [9, 10]. 


The potential flows of incompressible fluids with both / 1 5 

free (p = const, p*=0, where p is static pressure at the 
boundary, and p is the density of the medium adjacent to the 
subject flow) and straight rigid boundaries which are studied in 
these references include: 

a) flows which are stable or neutral to minor perturbations 
over the entire range of change in wave number of the imposed 
perturbed flow , except one isolated wave number, for which the 

flow is unstable. A typical feature of the configuration of j 

i 

these flows is the lack of a stagnation point in the finite j 

region of the main flow; 

b) flow which is stable or neutral to minor perturbations. 

The stagnation point of this flow lies on a rigid boundary of 
finite length; 

c) flow which is unstable for minor perturbations over a wide 
range of wave number changes and has a stagnation point within 
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the finite region of the main flow; 


d) flow which is unstable for minor perturbations over a wide 
range of wave number changes and has a stagnation point on a 
rigid boundary of infinite length in the finite region of the 
main flow. 

Thus, for the potential flows of incompressible fluids with 
partially free and partially rigid boundaries reviewed in these 
articles, instability occurs if the stagnation point lies in the 
finite region of the main flow and not on a rigid boundary of 
finite length. On the other hand, flows with a stagnation point 
o* tside the finite region of the main flow are virtually stable. 

This result of the linear theory of hydrodynamic stability is 
important for analyzing the causes of instability of flows with 
separation areas, which has been observed in experiments with 
specific ranges of R and M numbers for specific boundary 
conditions within the separation area. Since nonviscous flow in 
the 3 1 ream attachment region beyond the separation area which we 
have discussed is close in configuration to that class of flows 
with partially free and partially rigid boundaries studied in the 
theory of hydrodynamic stability, the conditions given above 
which determine the stability or instability of potential flows 
must be considered when analyzing the flow we are considering. 

The assumed position of the stagnation point in the 
nonviscous attachment region behind the separation area, a finite 
distance from the separation area in the plane of symmetry (or on 
the rigid boundary of infinite length when attachment is to the 
flat surface of an extended body beyond the nonsynnet r i c 
separation area) creates a situation in which potential analogs 
of the subject flow are unstable. On the other hand, we must 
bear in mind that potential analogs of the subject flow are 
stable, if nonviscous vortical flow in the attachment region 
beyond the separation area does not contain the stagnation point 
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in the finite region, and the configuration of the flow with an 
infinitely long backflow tongue develops. 


If we use the hypothesis that the instability of a nonviscous 
flow in an attachment region (caused, as shown by the examples of 
flows studied by the linear theory of hydrodynamic stability, by 
the configuration of flows in the attachment region) is the main 
reason for the instability of the entire viscous flow with 
separation area, then the rule for determining stable steady 
flows with separation areas will be that velocity in the finite 
region of the nonviscous vortical attachment not be zero. 

The validity of this hypothesis can be determined (at least 
at present) only by using the determination rule to compare the 
ranges of existence of stable stalled flows with the ranges of 
existence of stable steady flows observed in experiments over a 
wide range of changes in R and M numbers. Let us make this 
comparison . 

1. Incompressible fluid. According to calculations in points 
1-4 for flow around a flat slab with a separation area, the 
region of steady flows which satisfy the selection rule is 
limited by an R number, of about 1.2 x 10 2 . This R number 
corresponds to a strong dissipator in the separation area. If 
the dissipntor’s efficiency decreases, the R number at which a 
steady flow with separation area is disrupted diminishes. 

We know that experimental R numbers at which a steady flow is 
disrupted (called critical R^ number in literature), measured in 
wind or hydrodynamic tunnels are a function of the ratio of 
transverse dimension d of a cylindrical body to height II of the 
cross section of the working part of the tunnel. 

Thus, according to [11], steady flow around a round-section 
cylinder with a separation area is preserved up to R^=30 when 
d/H=0.025 and to R =G2 when d/H=0.!0. Therefore, calculated R. 
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values determined for an unrestricted current must be compared 
with experimental data obtained at small d/ll values. 4 With 
regard for thi3 circumstance, it follows from experimental data 
for flow around a round-section cylinder with a weak dissipator 
(friction on the back side of the cylinder) that R is about 30. 
If a dividing slab whose length is comparable to the length of 
the separation area is placed within an area with a plane of flow 
symmetry, steady flow is extended to higher H numbers. 

In this way Acrivos and his associates [12] successfully 

increased R $ several times to about 170 when d/H=0.05 (recent 

experiments by Acrivos* with small d/H values give an R^ of about 

$ 

100). Calculation of the force of friction applied to a dividing 
slab on the backward jet side shows that tractive force exceeds 
half the momentum of a backward jet whose parameters were 
calculated for a degenerate flow. 

A dividing slab is thus a strong dissipator, and the R^ of 
about 100 obtained by Acrivos* experiments is comparable to the 
predicted of about 120 for flow around n slab with a maximally 
strong dissipator. Base drag coefficients when R is 
approximately equal to are also comparable: calculated p^= 

-0.52, experiment p^=-0.45. 

Of course, when making this comparison we must keep in mind 
that the bodies differ in shape. However* the correlation 
obtained is quite promising for a first attempt, because, as far 
as we know, it has never been possible to obtain theoretical 
evaluations of the R^ numbers at which Karmen* s vortex street 
occurs since this flow chape was first observed. 


Calculated R ^ values for flow around a body located between 
parallel walls increase as d/ll increases. 

This comes to the author from V. V. Sychev. 
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The configuration of cavitation flows is close to that of the 
subject stalled flows. The fact that cavitation flows with 
closed cavities are steady flows at positive cavitation numbers Q 
=( p .-P k )/i/2pu^ is now commonly known [13]. However, using the 
determination rule for cavitation flows makes it possible to 
predict the existence of abroad region of stable steady 
cavitation flows when gas is released into a cavity and at 
negative cavitation numbers. In this case, given a nonviscous 
connection with various Bernoulli constants, there is generally 
no stagnation point beyond the cavity, or, given a negative 
cavitation number, Pq>P <<( on all streamlines coning from the 
cavity. 


2. Supersonic, velocities. Point 5 showed that, when R-~ and 
there is a strong .dissipator in the separation area, there are 
many possible steady flows which satisfy equations for gas 
dynamics of a viscous fluid, for which there are various 
corresponding base pressure values p^p^. However, only one, 
corresponding to point C (cf. figure 11) can be considered a 
stable steady flow, since the configuration of the flow in the 
nonviscous vortical attachment region which corresponds to point 
C satisfies the selection rule. 


As wo know, Chapman and associates [Id], in careful £1 

experiments on the stalled flow behind the reverse ledge on the 
surface of a body at 1.3<M<2.0 for R nuraberc from 0.6 x 10 4 to 
1.4 x 10 6 , observed the existence of laminar flow within a 
separation area at these high R numbers, in which base pressure 
varies as a function of fl number in quite close quantitative 
correspondence with values calculated for point C (p^ /p "1). 
Just as in Acrivos’s experiments, the surface of the body behind 
the ledge adjacent to the separation area is a strong dissipator 
in Chapman’s experiment, which justifies comparing experimental 
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data with calculations for a maximally strong dicsipator. 6 In 
summary, we can conclude that the hypothesis that configuration 
of flow in a nonviscous vortical attachment region has a decisive 
effect on the stability of steady flows with separation areas is 
confirmed by existing experimental data on laminar separation 
areas over a wide range of changes in R and M numbers. This 
makes it possible to apply this hypothesis in theoretical 
analysis of stalled flows. 

7. On the use of the selection rule in studying stalled 
flows . 

1* Turbulent low in a separation zone. As we know, the 
openness of the system of Reynold’s equations which describe 
vortical notion does not permit theoretical study of even the 
simplest instances of turbulent flows v/ithout using empirical 
quantitative data. In addition, use of empirical data makes it 
possible to calculate parameters for averaged turbulent flow for 
simple flow configurations (a wall turbulent boundary layer, a 
turbulent nixing layer on the boundary of a jet flowing out of an 
immersed space, the main portion of a turbulent jet). 


Note that Chapman’s use of the correlation between p n /p end 

up 

bottom pressure, observed in experiments at supersonic 
velocities (when flow was steady) and used in the Chaprzan- 
Korst calculating procedure (this procedure will be discussed 
in more detail later) for interpreting Roshko’s experimental 
data [15] on base pressure behind a round-section cylinder 
with n dividing slab at an R of about 10* 3 , is, in the light 
of what has been presented above, unjustified, since, in 
Roshko’s experiments, flow behind a round-section cylinder 
with a dividing slab was periodic (without constant 
properties) in nature at these R numbers. 
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We might expect that results satisfactory from an engineering 
standpoint can be obtained by calculation and by study of complex 
configurations of turbulent flow if they contain components which 
are simple configurations of turbulent flow studed in detail in 
"pure” form. 
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If these components cannot be identified in pure form in a 
turbulent flow of a complex configuration, then quantitative 
empirical data for the subject complex conf igurat ion of turbulent 
flow must be accumulated before calculation. Turbulent flow in a 
separation area has a complex configuration, but there is still 
little quantitative empirical data for it. Therefore, we will 
limit ourselves to discussing degenerate flow in which a 
turbulent mixing layer separates external flow from quiescent 
fluid within a separation area with constant pressure, since, in 
considering turbulent flow with a circulation core, we already 
cannot ignore turbulent friction stresses and energy dissipation 
within the core, which cannot as yet be calculated because of 
lack of quantitative empirical data for a random circulation core 
contour chape. 


For a completely turbulent flow in a separation area, R 
number determines only the nature of flow in the boundary layer 
on the surface of the body up to the flow separation point and 
the thickness of the boundary layer at the beginning of the 
turbulent nixing layer, which in turn define the position of flow 
separation on the body’s surface and indirectly affect U and 
values. These same equations, which describe flow in a turbulent 
nixing layer, do not contain molecular viscosity y and, 
consequently, R number. 


Let us consider flow around a slab at a 
equal zero and an M<<1 with degenerate flow within 
area from a maximally strong dissipator. 


which does not 
the separation 
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The decisive parameters of a turbulent nixing layer (U, /18 

Tp) are a function of ( 6 **) and empirical constant ct 
( proportional to the angle of expansion of the turbulent nixing 
layer). By making judgments similar to those in point 1 to 

Figure 13 


plotting c^=f(R) for a laminar nixing layer, we can, for a given 

(6%*), find possible steady flows for a turbulent degenerate flow 

which lie on curve BC on plane c^ , 1 /a (cf. figure 13). Point B, i 

as before, corresponds to limit constant-property flow (Pqp^B = ^‘ ' ' 

Obviously, only one point on this curve can represent real i 

flow, since a , the experimental constant for a turbulent j 

exp I 

nixing layer, is not a function of the shape of the body and ! 

separation area. j 

i ‘ 

i 

There are three possible cases. [ -j 

i - 

Case 1: l/o^. < ( 1/a ) ^ — r the assumed steady flow with : 

quiescent fluid within the separation area does not exist. Flow "■*: 

with a circulation core nay exist, but it cannot be calculated 
without regard for turbulent friction stresses and energy j 

dissipation within the core. 

i 

Case 2: l/o^^p) ( 1/a ) ^ — the assumed steady flow with 1 

quiescent fluid within the separation area does not satisfy the 

determination rule and cannot in fact be realized as steady. ? 

! 

Case 3: */° e xp = ^^ CT ^B — the t eady flow with 

quiescent fluid within the separation area exists and, since it 
satisfies the determination rule, can be realized in reality. 
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Since (l/o)g is a function of body shape, and l/° exp i- s 
constant, the probability that the proposed steady flow pattern 
with separation area will be observed in experiments with a 
random body shape is quite low. 


Thus, calculation of a turbulent separation area using a pure 
element — a turbulent mixing layer — is correct in those 
particular cases when the two conditions below are simultaneously 
satisfied: 

V° 7> < 1 /S) <f 1 /S e*P <7 - 1) 

If condition (7.1) is not fulfilled, flow in a separation 
area will be unsteady (i.e. it will possess increased pulsation 
amplitude in the turbulent nixing layer and a frequency 
characteristic different from that of an isolated (pure) 
turbulent mixing layer) or will remain steady, but will have a 
circulation core. In either case, calculation based on the 
concept of a pure nixing layer is unsuitable. 

In view of what has been presented here, it is obvious that 
the calculating procedure proposed by Korst [16] to describe 
actual turbulent flows with separation areas is incoiiect. Tnis 
procedure uses the concept of a pure turbulent mixing layer, but, 
of the three hydrodynamic equations which describe flow, only one 
— the equation of preservation of mass — is satisfied. 


Therefore, it is still absolutely necessary to satisfy 
conditions of type (7.1) to realize the proposed pattern of flow 
resulting from this complete system of hydrodynamic equations. 

It is conceivable that the well-known discrepancy between 
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That is, the dissipator’s characteristic is close to that of 
an ideal dissipator. 
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calculated data obtained using Horst’s method and data from 

systematic French [17] and English [18] experiments results 

because the flow pattern assumed in the calculation procedure 

cannot be realized in reality due to the fact that conditions of 

type (7.1) cannot be fulfilled. Note also that the failure /19 

to satisfy the momentum equation in Horst’s procedure leads to a 

paradoxical, but clearly unjustified, conclusion that we do not 

need to know of the value of the experimental constant o to 

exp 

calculate base pressure if there is turbulent flow in a 
separation area and no mass exchange between the separation area 
and the body. 

2 . Laminar flow in a separation area with unfixed flow 
separation point. Experiments on this important class show that 
flow in a laminar separation area is steady and has a circulation 
core over a wide range of R numbers up co about 10 6 if M>1. 

This flow has two regions for which calculations are done 
entirely differently: region I before the separation area to the 

flow separation point can be calculated (as is done in many 
studies using variations of integral ratio method) regardless of 
pressure within the separation area or taper; region II between 
the flow separation point and attachment, containing a 
circulation core, whose parameters are a function of the 
intensity of the descending jump or of taper. Calculating 
viscous flow in region II at high R numbers is a difficult 
problem whose solution requires effective methods which account 
for all the features of complex flow with sufficient accuracy, 
but which have not yet been found. We might expect that, as 
difficulties in numerical calculation of steady flow are 
overcome, we will be confronted by the problem of determining a 
single, practicable solution which must be reliable. 

In view of what has been presented here, it is conceivable 
that determination should also be made in terms of the 
configuration of flow in the nonviscous vortical attachment 
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region. To adapt the determination rule in this case, flow in a 

certain finite attachment region must be calculated as vortical 
and nonviscous. 


REFERENCES 


1. Lessen, M. , On the stability of the free laminar boundary 
layer between parallel streams, NACA, Rep. 979, 1949. 

2. Lock, R. C. The velocity distribution in the laminar boundary 
layer between parallel streams, Quart. J. Mech. and Appl. 
Math., 1951, vol. 4, No. 1, p. 42 

3. Efros, D. A. Hydrodynamic theory of flat parallel cavitation 
flows, Dokl . AN SSSR, 1946, vol. 51, No. 4. 

4. Gurevich, M. I., Teoriya struy ideal’noy zhidkc^ti [A Theory 
of Jets of an Ideal Fluid], Moscow, Fiznatgiz, 1961. 

5. Batchelor, G. E. , On steady laminar flow with closed 

streamlines at large Reynolds number, J. Fluid Mech. 1956, 
vol . 1 , pt . 2 . 

6. Batchelor, G. K. , A proposal concerning laminar wakes behind 
bluff bodies at large R numbers, J. Fluid Mech., 1956, vol. 

4 , pt . 4 , p . 386 . 

7. Lavrent ’yev, M. A., Var iat s ionnyy mctod v krayevykh zadachakh 
dlya s is ten uravneniy ell ipt icheskogo tipa [Variational 
Method in Extremum Problems for Systems of Elliptical Type 
Equations], Moscow, Academy of Sciences of the USSR Press, 
19G2 . 

8. Nikol’skiy, A. A., On flat vortical movements of gases, in 
IJekotoryye tochnyye resheniya uravneniy pros trans t vennykh 
techeniy gaca [Some Precise Solution to Equations of Three- 
dimenionsal Gas Flow], Tr. TsAGI, 1949. 

9. Ablow, C. M. , Hayes, V/. D., Perturbation of free surface 
flows. Technical Report No. 1, Office of Naval Research, No. 
7 onr-35807, Graduate Division of Applied Mathematics, Brown 
University, 1951 . 

10. Fox, J. L. , Morgan, G. M. , On the stability of some flows of 
an ideal fluid with free surfaces, Quart. Appl. Math. 1954, 
vol. 11, No. 4, pp. 439-456. 

11. Thoia, A., Flow post circular cylinders at low speeds, Proc. 
Roy. Soc. , 1933, A 141. 

12. Grove, A. S., Shair, F. H. , Petersen, E. E., Acrivos, A., An 
experimental investigation of steady separated flow past a 
circular cylinder, J. Fluid Mech., 1964, vol. 19, pt. 1, p. 
60. 



13. Sedov, L. I., Ploskiye zadachi gidrodinamiki i aerodinamiki 
[Flat problems in hydrodynamics and aerodynamics], Ed. 2, 
Moscow, "Nauka," 1966. 

D. R. , Kuehn, D. M. , Larson, H. K. , Investigation of 
separated flows in supersonic and subsonic streams with 
emphasis on the effe ct of transiti on. NACA Rep. 1356, 1958. 

15. Roshko, A., On the wake and drag of bluff bodies, J. 

Aeronaut. Sci., 1955, vol. 22, No. 2, p. 124. 

16. Korst, II. H. , A theory for base pressures in transonic and 
supersonic flow, J. Appl. Mech. , 1956, vol. 23, No. 4, pp. 
593-600. 

17. Carriere, P. , Recent studies by ONERA on problems of 
attachment, 7th Fluid Dynamic Symposium, Jurata, Poland, 

1965. 

18. Nash, J. F. , An analysis of two-dimensional turbulent base 
flow including the effect of the approaching boundary layer, 
ARC RM, 1963, No. 3344. 


39 




